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ABSTRACT 

Cygnus X— 2 appears to be the descendant of an intermediate-mass X-ray binary 
(IMXB). Using Mazzitelli's (1989) stellar code we compute detailed evolutionary se- 
quences for the system and find that its prehistory is sensitive to stellar input param- 
eters, in particular the amount of core overshooting during the main-sequence phase. 
With standard assumptions for convective overshooting a case B mass transfer starting 
with a 3.5M Q donor star is the most likely evolutionary solution for Cygnus X-2. This 
makes the currently observed state rather short-lived, of order 3 Myr, and requires a 
formation rate > 10 -7 — 10~ 6 yr _1 of such systems in the Galaxy. Our calculations 
show that neutron star IMXBs with initially more massive donors (ij 4M Q ) encounter 
a delayed dynamical instability; they are unlikely to survive this rapid mass transfer 
phase. We determine limits for the age and initial parameters of Cygnus X-2 and cal- 
culate possible dynamical orbits of the system in a realistic Galactic potential, given 
its observed radial velocity. We find trajectories which are consistent with a progenitor 
binary on a circular orbit in the Galactic plane inside the solar circle that received a 
kick velocity < 200 km s _1 at the birth of the neutron star. The simulations suggests 
that about 7% of IMXBs receiving an arbitrary kick velocity from a standard kick 
velocity spectrum would end up in an orbit similar to Cygnus X-2, while about 10% 
of them reach yet larger Galactocentric distances. 

Key words: accretion, accretion discs — binaries: close stars: evolution — stars: 
individual (Cygnus X-2) — X-rays: stars. 



1 INTRODUCTION 

The 9.84 d period neutron-star binary Cygnus X-2 has long 
been regarded as an archetypal long-period low-mass X- 
ray binary (Cowley et al. 1979, Webbink et al. 1983) where 
nuclear expansion of a Hayashi line giant drives mass trans- 
fer. Yet recent optical photometry and high-resolution spec- 
troscopy, while confirming that the donor has a low mass 
(Casares et al. 1998; Orosz & Kuulkers 1999), unambigu- 
ously showed that the spectral type of the optical counter- 
part is A5 ± 2 (Casares et al. 1998), significantly too hot for 
a Hayashi line donor. 

King & Ritter (1999; hereafter KR) and Podsiadlowski 
& Rappaport (2000; PR) argued that the system must be the 
descendent of an intermediate-mass X-ray binary (IMXB), 
with 3 — 4Mq as the likely initial donor mass. Such systems 
undergo a rapid mass transfer phase, previously regarded 
as fatal, with transfer rates exceeding the Eddington value 
by several orders of magnitude. For a narrow range of ini- 
tial separations the systems never reach the Hayashi line, 
or evolve away from it, during the subsequent phase with 
slower mass transfer. 



KR suggested that Cygnus X-2 is the product of a "case 
B" mass transfer sequence, where the donor star was already 
expanding towards the giant branch when mass transfer be- 
gan (Kippenhahn & Weigert [1967]). By contrast, PR pre- 
ferred an evolution where mass transfer started while the 
donor was still on the main sequence, with core hydro- 
gen burning terminating during the mass transfer phase. 
KR's considerations are semi-analytical and entirely based 
on generalised main-sequences (Giannone et al. 1968), while 
PR performed detailed binary sequences with a full stellar 
code. In Sect. 2 and 3 we use our evolutionary code to re- 
examine critically the case B solution rejected by PR, and 
resolve the discrepancy between KR and PR. To see how 
common a Cygnus X-2-like evolution is we consider IMXBs 
with still higher initial donor masses in Sect. 4. 

A further peculiarity of Cygnus X-2, addressed in 
Sect. 5, is its dynamical state. The system has a measured 
line-of-sight velocity of -208.6 km/s (Casares et al. 1998) 
and a Galactic latitude of -11.32° and longitude of 87.33°. 
At a distance of 11.6 kpc from the sun (Smale 1998), this 
places Cygnus X-2 at a Galactocentric distance of 14.2 kpc, 
and a distance from the Galactic plane of 2.28 kpc. Integrat- 
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ing the equations of motion in the Galactic potential we use 
Monte Carlo techniques to investigate possible trajectories 
of Cygnus X-2. 



2 CONSTRAINING THE PREHISTORY OF 
CYGNUS X-2 

The observed location of the donor star of Cygnus X-2 in 
the middle of the Hertzsprung gap in the HR diagram, and 
its small mass (0.4 — 0.7M©), suggest that the mass of the 
hydrogen-rich envelope remaining above the donor's helium 
core is very small (Mh/Mhc ^ 0.05). Stars with an initial 
mass 3— 4.5M0 leave the core hydrogen burning phase with a 
helium core mass Mh c in this range (e.g. Bressan et al. 1993). 
Mass transfer from such an intermediate-mass star on to a 
less massive neutron star is thermally unstable and involves 
an initial phase of rapid mass transfer with a highly super- 
Eddington transfer rate. This rapid phase lasts roughly un- 
til the mass ratio is reversed and the donor's Roche lobe 
expands upon further mass transfer. The subsequent slower 
transfer phase proceeds on the donor's thermal timescale, 
which is essentially given by the Kelvin-Helmholtz time of 
the donor when it left the main sequence. 

PR identified two possible routes from the post- 
supernova binary, i.e. after the formation of the neutron 
star, to the present system configuration. One route is via 
a genuine case B evolution, the other via an evolution they 
labelled 'case AB'. In the former the donor star has already 
left the main sequence when it fills its Roche lobe for the 
first time. The mass transfer phase is short-lived, with most 
of the mass transferred at a highly super-Eddington rate. 
Hence the neutron star mass increase is negligible, even if it 
accretes at the Eddington rate during the whole evolution. 
In the case AB solution discussed by PR mass transfer al- 
ready starts during the donor's main-sequence phase. The 
system briefly detaches when core hydrogen burning termi- 
nates, but mass transfer resumes when shell-burning is well 
established. As the rapid mass transfer phase terminates be- 
fore the donor leaves the main sequence the subsequent slow 
phase proceeds on a much longer timescale than in the gen- 
uine case B solution — the thermal time of the now less 
massive donor at the terminal main-sequence. The transfer 
rate is sub-Eddington for some time and the neutron star 
grows in mass. 

Given this we explore three different prescriptions to 
project the evolution of Cygnus X-2 backwards in time: 

(Model 1) To represent a case B evolution we assume 
that Mi =const., and that any material lost from the system 
carries the specific orbital angular momentum of the neutron 
star. Then the initial orbital separation en is 

M1+M2 ( M 2 \ 2 (2{M 2l -M 2 )\ 
a ^ a MTTWAM^) 6X H Mi J' W 

where M 2 i, M 2 denote the initial and present donor mass, 
and a the present separation (e.g. KR). This solution is valid 
only if a, > axivis, the orbital separation of a binary with 
a Roche-lobe filling terminal main-sequence (TMS) star of 
mass M 2 ,i. 

(Model 2) The same applies to a case AB evolution, 
except that we require en < otms- 

(Model 3) To account for a possible prolonged phase 



with sub-Eddington mass transfer in a case AB sequence 
we allow mass transfer to be conservative (total binary mass 
and orbital angular momentum is constant) for the last part 
of the evolution. We assume that the neutron star mass at 
birth was 1.4 Mq, and that the present neutron star mass is 
> 1.4 Mq. Hence the evolution backwards in time consists of 
two branches. During the conservative phase, characterised 
by aM 2 M 2 =constant, Mi reduces from its present value 
to 1.4 Mq. The second phase is calculated with constant 
neutron star mass (1.4M©), as in models 1 and 2 above. 

The present system parameters of Cygnus X-2 are: or- 
bital period Pq = 9.844 d, mass ratio q = 0.34±0.04 (Casares 
et al. 1998), and neutron star mass Mi = 1.4 — I.SMq. 
This mass range accommodates the canonical neutron star 
mass at birth, as well as the estimate by Orosz & Kuulkers 
(1999) based on modelling the observed ellipsoidal varia- 
tions. If we adopt a specific value Mi for the present neu- 
tron star mass, then the present donor mass is constrained 
to the range 0.30Mi M 2 ^ 0.38A/i by the observed mass 
ratio. This translates into a range of allowed initial separa- 
tions, as shown by the hatched regions in Fig. [l], for case 
B (model 1, wide spacing) and case AB (model 2, narrow 
spacing) evolution. The parameter space for case AB solu- 
tions with model 3 assumptions is tiny and hence not shown 
in the figure. For conservative mass transfer (model 3) the 
orbital separation of Cygnus X-2 today increases less steeply 
with decreasing donor mass than for the isotropic wind case 
(model 2). Hence the separation for model 3 sequences de- 
scribing the past evolution of Cygnus X-2 is generally larger 
than for the corresponding model 2 sequences. In particular, 
model 3 solutions require larger initial separations. Most of 
them are in conflict with the limit a; < otms, i.e. inconsis- 
tent with the assumption that the donor was on the main 
sequence when mass transfer started. 

More generally, if the mass lost from the system in se- 
quences with Mi = const, carries more (less) specific angular 
momentum than that of the neutron star, the required ini- 
tial separation is larger (smaller) than estimated by (1). A 
somewhat higher loss seems more likely than a smaller loss, 
hence this reduces the parameter space available for case AB 
solutions. 

Even if Fig. ^ indicates a viable solution it is not clear if 
the corresponding evolutionary sequence reproduces Cygnus 
X-2. Although the donor will have the observed mass (and 
hence radius) at the observed orbital period Po, the effective 
temperature is of course not constrained by the above con- 
siderations. To check this we need a detailed binary sequence 
with full stellar models, as presented in the next section. In 
particular, it appears that the case AB sequences described 
by PR require the donor star to be already fairly close to the 
end of core hydrogen burning at turn-on of mass transfer. 
This is likely to narrow the parameter space for case AB 
sequences even further. 

We note that the limits shown in Fig. |l] depend via htms 
somewhat on the stellar input physics, in particular on the 
amount of convective core overshooting (cf. the discussion 
in the next section). 
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Figure 1. Allowed range for the initial orbital separation a; 
as a function of initial donor mass M-ii, for three different fixed 
neutron star masses, as labelled. The hatched areas correspond 
to case B evolution (model f ; wide spacing) and case AB evolu- 
tion (model 2; narrow spacing), assuming constant neutron star 
mass. The parameter space for case AB solutions with model 3 
assumptions is so small that it is not shown. The almost hori- 
zontal boundary at Oi = citms — 11 -Reg separating case B and 
case AB solutions represents the terminal main sequence, iitms 
was calculated using main— sequence models by Dominguez et al. 
(1999). See text for details. 

3 MODEL CALCULATIONS 

We calculated several detailed early massive case B binary 
evolution sequences, using Mazzitelli's stellar code in its 
1989 version (see Mazzitelli 1989, and references therein, for 
a summary of the input physics) . Mass transfer was treated 
as in Kolb (1998). Table [l] summarises initial and final sys- 
tem parameters. The neutron star mass was 1.4M© = const, 
in each case, with angular momentum loss treated as in 
model 1 above. Specifically, we present in detail the evo- 
lution of a system with initial donor mass M2 = 3.5M0 and 
orbital separation at = 13.3Rq (Sequence I) and 11.17?© 
(sequence II). 

The figures confirm the general behaviour of case B 
evolutionary sequences as described above. Sequence I gets 
very close to the observed state of Cygnus X-2, although it 
fails to fit all observed parameters simultaneously with high 
accuracy. The evolutionary track in the HR diagram passes 
through the Cygnus X-2 error box, but at a period slightly 
longer than the observed value Po = 9.844 d. This can be 
seen in Fig. EL which shows that at Po the sequence I donor 



is somewhat too cool, while in sequence II it is somewhat too 
hot. Hence the initial separation for a simultaneous fit of P, 
T e ff (and L) lies between 11.17?© and 13.37?0. The transfer 
rate at Po for sequence I is 5 x 10~ 8 M©yr~ 1 , a factor 2 — 3 
higher than the Eddington rate. 

Observational estimates place the actual accretion rate 
in Cygnus X-2 close to the Eddington limit. With no ev- 
idence of significant mass loss from the system we expect 
that the transfer rate is also of the order of the Eddington 
rate. As the transfer rate in the model sequence is set by 
the thermal time of the progenitor star, a somewhat smaller 
initial secondary mass (e.g. 3.OM0) would give a lower rate; 
this has already been noted by PR. 

We note the following main differences between our se- 
quence I and the case B sequence calculated by PR: 

(1) PR chose 21.37?© as the initial separation. This leads to 
the combination 1.4 + O.66M0 for the component masses at 
P = Po, i.e. to a mass ratio 0.47, significantly larger than 
observed. Hence the PR case B sequence represents a poor 
fit for Cygnus X-2. 

(2) More importantly, the mass transfer rate in the 
slow phase of sequence I decreases from 10 _7 M©yr _1 to 
10 _8 M©yr _1 over a period of 2.5 Myr, while in PR's model 
the transfer rate is always larger than 10 M©yr _1 (which 
is inconsistent with observations). A similar calculation by 
Tauris et al. (2000) with essentially the same stellar code as 
the one used by PR gives a similarly high rate in the slow 
phase. 

The initial separation for our sequence II is essentially 
the same as the one used by PR for their case AB evo- 
lution (11.57?©, corresponding to a stellar radius 5.6P© of 
the donor at turn-on of mass transfer). The most impor- 
tant difference between the two calculations is the degree of 
convective overshooting assumed during the main-sequence 
phase: none in our models, a very strong one in PR's models. 
Hence PR's main-sequence band is much wider, with 6.6P0 
as the maximum radius J?2,max that a 3.5M0 star reaches 
during core hydrogen burning, compared to 3.67?© in our 
case. A moderate extent of core overshooting is favoured in 
the literature, giving J?2,max = 5.07?© (Schaller et al. 1992), 
5.77?© (Bressan et al. 1993) and 4.87?© (Dominguez et al. 
1999) for a 3.5M0 star with solar composition. Schroder 
et al. (1997) — who use the same code and input physics 
as PR — prefer a rather efficient overshooting leading to 
values up to 7?2,max = 7.17?©, while in a subsequent pa- 
per (Pols et al. 1997) the same authors concluded that 
7?2,max = 5.97?© (which corresponds to their overshooting 
parameter 8 — 0.12) provides the best overall fit to observa- 
tions. This is only barely larger than the donor's radius at 
turn-on of mass transfer in PR's case AB solution. Hence 
standard assumptions about overshooting favour a case B 
solution for Cygnus X-2 over a case AB solution. 

We conclude that a case B solution is a viable fit for the 
evolutionary history of Cygnus X-2. Given the parameter 
space limitations for a case AB evolution it does appear as 
the more likely solution for Cygnus X-2, even though this 
implies that the presently observed state of Cygnus X-2 is 
rather short-lived, of order several million years. 

PR pointed out that the surface composition of the 
donor in Cygnus X-2 would be significantly hydrogen- 
depleted (X ~ 0.1) and show signs of CNO-processing if 
it had undergone a case AB evolution. Unfortunately, this 
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Table 1. Donor mass M 2 f and orbital period Pf when mass 
transfer ceases for sequences with different initial donor mass 
M 2 i and initial orbital period P; (initial orbital separation a,). 
Sequences labelled "runaway" encounter the delayed dynamical 
instability; M 2 / and Pf refer to the system at this point. 



2.5 





1 1 I 1 1 1 


i 

/ - 










\ \ - 5 






M 2 /M e = 3.o\ \ 


0.6 


, i , 


2.o\ 

, , 1 , , , 1 , 


1.0 - 



4.2 4 3.8 3.6 



log T e „ [K] 



Figure 2. HR diagram showing sequence I (solid) and a track of 
a 3.5M0 single star (dotted), calculated without core overshoot- 
ing. Labels along the track indicate the donor mass. The shaded 
error box marks the observed location of Cygnus X-2 (lower and 
upper bound correspond to donor masses 0.4Mq and O.7M0, re- 
spectively) . 



does not distinguish unambiguously between a case AB and 
case B evolution, as the same is in principle true for a donor 
that had undergone a case B mass transfer. The models in 
our sequence I close to the position of Cygnus X-2 have a 
surface hydogen mass fraction of X = 0.29, while C/N and 
O/N are close to the equilibrium values for CNO burning at 
2 x 10 7 K. 
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Figure 3. Mass transfer rate M vs time (elapsed since the be- 
ginning of core hydrogen burning) for sequence I. Wiggles reflect 
numerical noise due to the explicit treatment of mass transfer. 
The instabilities arc particularly strong at the end of the rapid 
phase close to t = 173.8 Myr. The semi— regular dips of M during 
the slow phase are due to sudden changes of the stellar radius, 
caused by a discontinuous change of the inner boundary of the 
convective envelope. 



4 SEQUENCES WITH INITIALLY MORE 
MASSIVE DONOR STARS 

The calculations discussed so far represent early case B mass 
transfer solutions for systems with a moderate initial mass 
ratio, q = M 2 /Mi <; 2.5. Systems with smaller mass ratio 
(q ^ 1) have no initial rapid mass transfer phase and can 
be understood semi-analytically (see e.g. Kolb 1998, Ritter 
1999), while the fate of systems with yet larger initial mass 
ratio <; 3 is unclear. Hjellming (1989) pointed out that sus- 
tained thermal-timescale mass transfer can lead to a delayed 
transition to dynamical-timescale mass transfer. The reason 
for this is that the adiabatic mass-radius index of initially 
radiative stars decreases significantly when the radiative en- 
velope is stripped rapidly. This causes the Roche lobe to 
shrink faster than the star. Estimates for the critical initial 
mass ratio which just avoids the delayed dynamical insta- 
bility give q a 3 — 4 (Hjellming 1989; Kalogera & Webbink 
1996), although none of these are based on self-consistent 
mass transfer calculations. 

This critical value is relevant for identifying possible de- 
scendants of neutron-star systems undergoing early massive 
case B mass transfer. These should appear as binary mil- 
lisecond pulsars lying significantly below the orbital period- 
white dwarf mass relation found by Rappaport et al. (1995; 
see also Tauris & Savonije 1999 for calculations with updated 
input physics) for systems descending from Hayashi line 
low-mass X-ray binaries. Obvious candidates are systems 
with a fairly high-mass white dwarf but short orbital pe- 
riod. The most discrepant system is B0655+64 (P = 1.03 d, 



© 0000 RAS, MNRAS 000, 000-000 



The violent past of Cygnus X-2 5 




4 6 8 10 12 
P [d] 



Figure 4. Upper panel: Mass transfer rate M vs orbital period P 
for sequence I. The numerical noise at the end of the rapid phase 
(7 < P/d < 8) is suppressed. The dotted vertical line marks the 
period of Cygnus X-2. Lower panel: Effective temperature vs P, 
for sequence I (solid) and II (dashed). The filled and open circle 
on the tracks mark where the donor mass is 0.6Mq and 0.5Mq, 
respectively. The short bar indicates the location of Cygnus X— 2. 



Mwd <; 0.7; cf. KR). The fairly massive white dwarf implies 
a donor mass 5M0 in the progenitor binary. 

We performed test calculations similar to our sequences 
I and II, but with initial donor mass 3.75Mq, 4.OM0 and 
5.OM0. While the 3.75Mq sequence was stable through- 
out, the 4.OM0 and 5.OM0 sequences encountered runaway 
mass transfer (where we stopped the calculations) rather 
early in the rapid mass transfer phase (see Tab|l|). An ad- 
ditional test sequence with a 3M© donor near the end of 
core hydrogen burning and initial mass ratio 4 (i.e. pri- 
mary mass O.75M0, assumed constant, as above) encoun- 
tered runaway mass transfer at a donor mass 2.6M0. This 
is in perfect agreement with Hjellming's prediction, based 
on his Fig. IV. 1 and a Roche lobe curve corresponding to 
g = 4. 

Unless a neutron star binary can survive even such a 
dynamical-timescale mass transfer, the phase space avail- 
able for a Cygnus X-2-like evolution is severely limited by 
this upper limit on the initial donor mass. Our calcula- 
tions suggest that if most neutron stars form with 1.4Mq 
the maximum white dwarf mass in an endproduct of early 
massive case B evolution is O.55M0. If this is true, nei- 
ther B0655+64 nor the recently discovered system J 1453-58 



(P = 12.422 d, M WD ^ 0.88M Q ; cf. Manchester et al. 1999) 
can have formed in this way. 

The occurrence of the delayed dynamical instability is 
intimately linked to the internal structure of the donor star. 
Therefore it is not surprising that the maximum initial donor 
mass for early case B mass transfer, just avoiding this in- 
stability, depends on the stellar input physics. Using an up- 
dated version of Eggleton's stellar code (see e.g. Tauris & 
Savonije 1999), Tauris et al. (2000) found 5 M for this limit, 
with a corresponding maximum white dwarf mass 0.8 Mq in 
millisecond pulsar binaries formed in this way. 



5 THE TRAJECTORY OF CYGNUS X-2 

As described in the introduction, Cygnus X-2 is in an un- 
usual dynamical state. In this section we investigate possible 
trajectories of Cygnus X-2, given the observed line-of-sight 
velocity and the constraints on the evolutionary state as de- 
rived in Sees. 2 and 3. 

In Fig. 5 we plot the cumulative velocity distribution 
for systems with initial separations of d = 15Rq, i.e. imme- 
diately after circularisation of the binary orbit following the 
SN explosion producing the neutron star. The secondary 
is taken to have a mass M2 = 3.5 Mq while the mass of 
the helium star prior to the supernova is Mi = 5.0 Mq. 
We consider two different distributions for the kick velocity 
imparted to the neutron star; namely those by Hansen & 
Phinney (1997) and Fryer (1999). Throughout the following 
section, when we refer to kick velocities we mean the kick 
imparted on the binary allowing for the effects of mass-loss 
and an asymmetric supernova, and not simply the kick im- 
parted on the neutron star from the latter. 

Given its present position and velocity, we integrate the 
trajectory of Cygnus X-2 backwards to the birth of the neu- 
tron star using a model for the Galactic potential suggested 
by Paczyriski (1990) (reviewed in the Appendix). From sec- 
tions 2 and 3 it is clear that the age t of the neutron star is 
essentially the main-sequence lifetime of the progenitor star 
of the present donor. To account for the allowed range of the 
initial donor mass and for uncertainties from the width of 
the main-sequence we choose 150 < t/Myr < 275. As only 
the line-of-sight velocity, v; os , is known, we considered a set 
of trajectories where the current velocity was given by 

V = Vq + Vl os + QfVi + I3-V2 (2) 

where Vi os , vi, and V2 are mutually orthogonal, and vi has 
a zero component in the z direction. When computing tra- 
jectories, it is most convenient to work in units related to 
the mass and size of the galaxy. As a natural unit of velocity 
we use 207 km/s, the Kepler velocity at 1 kpc distance from 
a point mass 10 10 Mq. In the above equation vi, and V2 are 
given unit lengths in our code units. 

We integrated the equations (A9) for a range of values 
of a and (3 and noted the radius, r, and time, t, whenever 
the trajectory cut the Galactic plane. We also computed the 
kick velocity, Vk the system must have had if the neutron star 
has formed at that time, assuming the system had previously 
been in a circular orbit in the Galactic plane. 

In Fig. 6 we plot the values of r and Vk for all values 
of a and f5. We note that for r <J 10 kpc, 140 km/s <J Vk 
600 km/s. This range of required values of Vy should be 
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Figure 5. The cumulative velocity distribution (in km/s) for 
systems with initial separations, d = 15J?0. This is the velocity 
of the centre-of-mass of the binary as opposed to the kick of 
the neutron star itself. Dotted line is computed from the neutron 
star kick distribution of Hansen & Phinney (1997), whereas the 
dashed line is for the kick distribution of Fryer (1999). 



Figure 6. The kick velocity, Vk (in km/s), of binaries as a func- 
tion of formation radius, r (in kpc). In the case of the larger dots, 
150 < t < 275, and 5 < r < 10. 



compared with that expected given the neutron star kick 
distributions of Hansen & Phinney (1997) and Fryer (1999) 
in Fig. 5. It is clear from Fig. 5 that the expected value of Vk 
is far lower than that required for many of the trajectories 
shown in Fig. 6. In other words, in a large fraction of systems 
formed at a radius r ^ 10 kpc, the trajectory of the binary 
will confine the system to the inner regions of the Galaxy; 
Cygnus X-2 must have been relatively unusual in reaching 
out beyond the solar circle. 

We plot the time of formation, t, as a function of r for 
all values of a and (3 in Fig. 7. As in Fig. 6, the data points 
where 150 Myr < t < 275 Myr and V k < 200 km/s are plot- 
ted as larger dots. From the work of Sections 2 and 3 these 
are the values most likely to be applicable to a progenitor of 
Cygnus X-2. We note from Fig. 7 that the binary is unlikely 
to have originated within 5 kpc of the Galactic centre as tra- 
jectories cutting the Galactic plane at such small radii do 
so at the wrong time. Fig. 6 also shows they would require 
unreasonably large kicks. 

In Fig. 8 we plot the values of a and f3 for those tra- 
jectories which cut the Galactic plane t Myr ago, where 
150 < t < 275, with system kick velocities required to take 
the binary from a circular orbit restricted by Vk < 200 km/s. 
Values of a and (5 besides those shown in Fig. 8 were consid- 
ered, but none satisfied the above conditions. We see that 
only a small number of possible trajectories are credible with 
these restrictions applied. Here we note the effect of the 
Galactic rotation. For a <, 0.4 the trajectories are in the 
opposite direction to the Galactic rotation. In other words 
the kick the system receives on formation will oppose its ini- 
tial velocity in a circular orbit, and will therefore have to be 
larger than for a similar trajectory in the same direction as 
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Figure 7. The age of formation, t (in Myr), of binaries as a 
function of formation radius, r (in kpc). For the larger dots, 150 < 
t < 275, and V k < 200 km/s. 



the Galactic rotation. Here the effect of the kick is boosted 
by the initial velocity of the circular orbit. 

For illustration, in Fig. 9 we plot the trajectory (R and 
z) for a = 0.75 and f3 = —0.2. In this case, the trajectory 
cuts the Galactic plane on two occasions within the last 
250 Myr, the first occasion being 95 Myr ago, at a radius, 
Rcut = 14.1 kpc, and the second occasion being 172 Myr 
ago, at a radius , R cut = 5.9 kpc. If the system was formed 
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Figure 8. The values of a and j3 producing trajectories where 
the binary was ejected from the Galactic plane from within a 
radius, r < 10 kpc. The system is assumed to have received a 
kick, Vk < 200 km/s, and the ejection to have occurred a time t 
Myr ago, where 150 < t < 225 (open squares) and 225 < t < 275 
(filled squares). 




5 10 15 20 

R 

Figure 9. Side view of the Galaxy (its plane is the line z = 0), 
showing the trajectory of the binary for a = 0.75 and fi = —0.2. 
The dashed curve shows the path from formation to the present 
day, and the solid curve shows the path that the system will follow 
over the next one billion years. 



Figure 10. The Galactocentric distance (in kpc) reached by the 
binary as a function of time (in Myr) for 100 randomly-picked 
initial conditions. 

in the disc at this latter time, and was initially travelling on 
a circular orbit, it would have had to receive a kick, VL ~ 130 
km/s, which sits comfortably within the range given in Fig. 
5. 

Thus far we have demonstrated that there are viable 
trajectories for Cygnus X-2 compatible with the observed 
line-of-sight velocity. We now ask the reverse question: what 
fraction of binaries originating at some radius will produce 
systems resembling Cygnus X-2? Given that we expect sys- 
tems to have originated within the solar circle, as most mas- 
sive stars are located within 10 kpc of the Galactic centre, 
we might suspect that Cygnus X-2 is relatively unusual in 
being at a larger distance from the Galactic centre and sig- 
nificantly away from the Galactic plane. By Monte Carlo 
simulation we were able to produce a large number of sys- 
tems at various initial radii, assumed to be initially on circu- 
lar orbits, with kick velocities drawn from the distributions 
given in Fig. 5, and follow their trajectories in the Galactic 
potential. Investigation demonstrated that in determining 
whether a particular trajectory would produce a Cygnus X- 
2 like system, the maximum radius reached by the system, 
i?max was a good diagnostic. This is illustrated in Fig. 10 
where we plot the radii of 100 binaries as a function of time 
where the initial trajectories are drawn randomly, from the 
kick distribution of Hansen and Phinney (see Fig. 5), and 
choosing the initial radius to be between 2 kpc and 10 kpc. 
The trajectories fall into two categories: those where the bi- 
nary remains at a radius similar to that at which it was 
located initially, and those where the binary is ejected sig- 
nificantly into the Galactic halo, reaching maximum radii of 
<; 10 kpc. Cygnus X-2 clearly belongs to the latter category. 
In order for a system to be at a radius today similar to that 
of Cygnus X-2, we require 15 <J i? m ax ^ 20 kpc, providing 
the system originated in the Galactic disc somewhere within 
10 kpc of the Galactic centre. A number of systems will be 
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Figure 11. The maximum Galactocentric distance reached by 
the binary, i? m ax, as a function of formation radius, R. One hun- 
dred different trajectories are considered here. The initial radius 
was chosen randomly from two to ten, and the kick velocity drawn 
from the distribution given by Hansen and Phinney (see Fig. 5). 



Figure 12. The fraction of systems on similar trajectories to 
Cygnus X-2 (solid line), and those on longer period orbits (dashed 
line) as a function of formation radius, R (in pc). 



ejected to even larger radii. In such cases a binary of age 
similar to Cygnus X-2 would still be on the outward bound 
portion of its orbit. 

In Fig. 11 we plot i? max as a function of initial radius 
R. In this case the initial radius was chosen randomly from 
2 kpc to 10 kpc, and the kick velocity drawn from the dis- 
tribution given by Hansen and Phinney (see Fig. 5). From 
this figure we note that Cygnus X-2-like objects may be 
produced when r <; 5 kpc, and that the relative frequency is 
relatively independent of formation radius. Systems having 
larger values of 7? max will also be produced, the frequency 
increasing with radius of formation. The relative frequency 
for Cygnus X-2-like systems and those on longer period or- 
bits is plotted as a function of R in Fig. 12. 

Assuming that the formation rate of binaries scales 
as the surface density of stars, which is given by E ~ 
E e" r/3kpc (see Binney & Merrifield 1999), we computed the 
relative number of Cygnus X-2-like systems which would be 
produced by integrating over the entire disc within the so- 
lar circle. The fraction of systems resembling Cygnus X-2 
(/ cyg ) and the total fraction of systems which will travel sig- 
nificantly further out than their formation radius (/out) are 
listed in Table ^ as a function of the initial separation within 
the post-supernova binary, d, the primary helium star mass 
prior to the supernova, Mi , and the mass of the secondary, 
M2. We found that ~ 7% of all binaries will produce Cygnus 
X-2-like binaries on trajectories that would place them at 
a radius similar to Cygnus X-2 today. A further ~ 5 — 15% 
will be further out than Cygnus X-2, whilst the remainder 
will be located closer to the Galactic centre. These results 
apply equally to both velocity distributions plotted in Fig. 5. 



d/Re 


Mi/M@ 


M 2 /M Q 


/cyg 


/out 


10.0 


5.0 


3.5 


0.07 


0.17 


15.0 


5.0 


3.5 


0.06 


0.15 


20.0 


5.0 


3.5 


0.06 


0.13 


25.0 


5.0 


3.5 


0.06 


0.11 


15.0 


5.0 


4.5 


0.06 


0.13 


20.0 


5.0 


4.5 


0.06 


0.12 


30.0 
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4.5 


0.05 


0.09 


10.0 


4.0 


3.5 


0.06 


0.13 


10.0 


6.0 


3.5 


0.08 


0.23 


15.0 


4.0 


4.5 


0.05 


0.09 



Table 2. The fraction of systems resembling Cygnus X-2 (/cyg) 
and the fraction of systems which will be located further out 
than Cygnus X-2 (/out) as a function of the initial separation 
within the post-supernova binary, d, the primary mass prior to 
the supernova, M 1 , and the mass of the secondary, M2 . 

6 SUMMARY 

Our binary evolution calculations with full stellar models 
have verified that Cygnus X-2 can be understood as the 
descendant of an intermediate-mass X-ray binary (IMXB). 
The most likely evolutionary solution is an early massive 
case B sequence, starting from a donor with mass ~ 3.5M0. 
This implies that the presently observed state is rather 
short-lived, of order 3 Myr. This in turn points to a large 
IMXB formation rate of order 1(T 7 - 1CT 6 yr _1 if Cygnus 
X-2 is the only such system in the Galaxy. It is likely that 
there are more as yet unrecognized IMXBs, so the Galactic 
IMXB formation rate could be even higher. Hercules X-l 
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seems to be in the early stages of a similar case A or case 
AB evolution (e.g. van den Heuvel 1981). The prehistory of 
Cygnus X-2 is sensitive to the width of the main-sequence 
band in the HR diagram, i.e. to convective overshooting in 
that phase. The alternative evolutionary solution for Cygnus 
X-2 suggested by PR, a case A mass transfer followed by a 
case B phase, is viable only if overshooting is very effective, 
i.e. more effective than hitherto assumed in the literature. 

Using Mazzitelli's stellar code (Mazzitelli 1989) we 
found that neutron star IMXBs that start case B mass trans- 
fer with initial donor mass > 4M0 will encounter a delayed 
dynamical instability. (Note that the maximum donor mass 
that just avoids this instability depends on stellar input 
physics). The components are likely to merge and perhaps 
form a low-mass black hole. Given the high formation rate 
of Cygnus X-2-like objects the formation rate of such black 
holes could be rather substantial. 

We have shown that the large Galactocentric distance 
of Cygnus X-2 and its high negative radial velocity do not 
require unusual circumstances at birth of the neutron star. 
There are viable trajectories for Cygnus X-2 that firstly cut 
the Galactic plane inside the solar circle at a time consis- 
tent with the evolutionary age of Cygnus X-2, and secondly 
which require only a moderate neutron star kick velocity 
(< 200km/s). We estimate the fraction of systems resem- 
bling Cygnus X-2, i.e. of systems on orbits that reach Galac- 
tocentric distances > 15 kpc, as 10 — 20% of the entire IMXB 
population of systems. 
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APPENDIX: THE GALACTIC POTENTIAL 

The Galactic potential can be modelled as the sum 
of three potentials. The spheroid and disc components are 
given by 

GM a 



1>s{R,z) 



(R 2 + [a s + (z 2 + b 2 y/ 2 } 2 ) 1/2 

GMd 

(R 2 + [a d + (z 2 + b 2 y/ 2 } 2 ) 1/2 



(1) 
(2) 



where R 2 — x 2 + y 2 . The component from the Galactic 
halo can be derived assuming a halo density distribution, 
p h , given by 

Pc 



Ph = 



1 + (r/r c ) 2 



where r — x + y + z 
yields the potential 



(3) 

The above density distribution 



GM C 



5* (i + J) + 7 atan ik) 



(4) 



where M c = AiYp c r 3 . The total Galactic potential is the sum 



(3) 



Following Paczyhski (1990), we use the following choice of 
parameters: 

a s = 0, b s = 0.277 kpc, M s = 1.12 x 10 10 M Q , (6) 

a d = 3.7 kpc, b d = 0.20 kpc, M d = 8.07 x 1O 1O M , (7) 

r c = 6.0 kpc, M c = 5.0 x 10 10 M Q , (8) 

Because of the cylindrical symmetry of the potential, the 
integration of the trajectories can be simplified to consider 
the evolution of the z and R components only, as given below 



dR dz 

It = VR > dt =V * 



dVR 

dt 



KdRl, 



+ 



3z 



dv z 



R 3 ' dt 



\dz) R 



0) 



where the z component of the angular momentum, j z 
Rvtp. 
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